The dynamics of laser droplet generation 
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We propose an experimental setup allowing for the characterization of laser droplet generation in terms of 
the underlying dynamics, primarily showing that the latter is deterministically chaotic by means of nonlinear 
time series analysis methods. In particular, we use a laser pulse to melt the end of a properly fed vertically 
placed metal wire. Due to the interplay of surface tension, gravity force and light-metal interaction, undulating 
pendant droplets are formed at the molten end, which eventually completely detach from the wire as a conse- 
quence of their increasing mass. We capture the dynamics of this process by employing a high-speed infrared 
camera, thereby indirectly measuring the temperature of the wire end and the pendant droplets. The time series 
is subsequently generated as the mean value over the pixel intensity of every infrared snapshot. Finally, we 
employ methods of nonlinear time series analysis to reconstruct the phase space from the observed variable and 
test it against determinism and stationarity. After establishing that the observed laser droplet generation is a 
deterministic and dynamically stationary process, we calculate the spectra of Lyapunov exponents. We obtain a 
positive largest Lyapunov exponent and a negative divergence, i.e., sum of all the exponents, thus indicating that 
the observed dynamics is deterministically chaotic with an attractor as solution in the phase space. In addition 
to characterizing the dynamics of laser droplet generation, we outline industrial applications of the process and 
point out the significance of our findings for future attempts at mathematical modeling. 
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The dripping faucet is one of the paradigmatic examples 
of deterministic chaos. Due to the inherently nonlinear in- 
terplay between the surface tension, mass of the droplets 
and the dripping rate, the system exhibits extreme rich- 
ness of dynamics, culminating in the emergence of deter- 
ministically chaotic behavior. Here we propose and exam- 
ine the dynamics of a conceptually closely related process, 
which however, is governed by additional physical phe- 
nomena. While the complexity of laser droplet genera- 
tion also relies on the interplay between the surface ten- 
sion and the increasing droplet mass, the addition of light- 
metal interaction and the fact that molten metal has differ- 
ent properties than water warrant diversity if compared 
to the traditional dripping faucet experiment. It is there- 
fore all the more fascinating that from the viewpoint of dy- 
namics, the two processes share deterministic chaos as an 
inseparable ingredient. Indeed, our analysis reveals that 
the phase space, reconstructed from the indirect tempera- 
ture measurements of the metal droplets, is characterized 
by an attractor having negative divergence and a positive 
largest Lyapunov exponent. As such, it has all the prop- 
erties that are characteristic for deterministically chaotic 
systems. These essential insights into the dynamics of laser 
droplet generation from experimental data are paramount 
for the proper introduction of the process to actual indus- 
trial applications, as well as to modeling attempts that may 
further facilitate its understanding. Demonstrating the 
emergence of chaos in a realistic engineering setup adds to 
the evergreen nature of the subject, and in this sense, we 
hope that the study will be inspirational and spawn fur- 
ther research aimed at unraveling the dynamics of laser 
droplet generation and revealing its full potentials. 



I. INTRODUCTION 

Nonlinear dynamical systemsi offer a gateway to fascinat- 
ing phenomena that imbue many facets of our existence. Al- 
though frequently going by unnoticed, deterministic chaos, 2 
fractal structures, 3 synchronization, 4 and even the stochastic 
resonance, 5 are phenomena that are at the very heart of nu- 
merous manmade and natural systems. Be it only a thought 
or a heart beat, 6 the organization of traffic, 7 or the weather 
forecast, 8 nonlinear dynamics plays an important role in it all. 
However, while the dripping faucet 9 or the flapping of butter- 
fly wings 10 are paradigmatic examples of deterministic chaos, 
their omnipresence and universal appeal are in stark contrast 
with many of the processes in engineering and technical sci- 
ences, where the complexity underlying them frequently re- 
mains unexplored or at least unknown to the wider audience. 
Since mathematical models for complex processes are diffi- 
cult to construct, and are therefore either non-existent or cap- 
ture only the essential ingredients of the dynamics, one of the 
common obstacles to overcome is the characterization of the 
process from experimental or observed data sets. 

Nonlinear time series analysis 11 enables the determination 
of characteristic quantities, for example the number of ac- 
tive degrees of freedom or invariants such as the Lyapunov 
exponents, 2 of a particular system by analyzing the time 
course of one of its measured variables. In the past two 
decades numerous successful applications of nonlinear time 
series analysis on data sets from the most diverse fields of 
research have been reported, 12 and there are still new ap- 
proaches being proposed to date. One of the most recent ad- 
vances is the merging of concepts from the theory of com- 
plex networks 13 and time series analysis, 14 particularly also 
recurrence plots, 15 which together gave rise to new quantifiers 
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for experimental data sets. Despite the availability of com- 
prehensive and well-documented programs 17 , however, there 
are still branches of science in which the application of these 
methods could lead to substantial advancements. The failure 
of applying them is in part certainly due to the difficulties 
associated with successfully bringing together scientists that 
are working on very different and seemingly completely dis- 
joint subjects, but also due to the fact that not all data sets are 
equally amenable to methods of nonlinear time series analy- 
sis. While the latter indeed offer tools that bridge the gap be- 
tween experimentally observed irregular behavior and the the- 
ory of dynamical systems, it should be emphasized that this is 
true foremost if the series under study has properties that are 
typical of deterministic dynamical systems.— ~— Moreover, is 
has to be verified if the observed irregular behavior originated 
from a stationary system, 21 ' 22 for it may solely be a conse- 
quence of varying system parameters during data acquisition. 
These are important issues that have to be addressed, espe- 
cially on experimental recordings, as we will try to emphasize 
throughout this work. 

In this paper we propose an experiment that allows us to 
characterize the process of dripping via laser-induced heat- 
ing of the end of a properly fed metal wire, i.e., the laser 
droplet generation. Particularly, we are interested in the dy- 
namical properties of this process, which via the analogy with 
the dripping faucet 9 promise to be very interesting. However, 
although being conceptually similar, we note that the laser 
droplet generation is governed by additional physical phe- 
nomena. While the surface tension and droplet mass also play 
a crucial role, the additional effects brought about by light- 
metal interaction, heating, phase transitions and the fact that 
molten metal has different properties than water distinguish 
the process significantly from the traditional dripping faucet 
experiment. In our case a laser pulse is used to melt the end of 
a vertically placed metal wire. Due to the interplay between 
the surface tension and the gravity force a pendant droplet is 
formed from the molten end. The pendant droplet eventually 
becomes fully detached due to its growing mass and the laser 
light-matter interaction. During the process of droplet for- 
mation the wire has to be properly fed, thereby mimicking 
the flow rate of water in the dripping faucet experiment. The 
most important variable to monitor during the process is the 
temperature of the wire end and the pendant droplet, which 
we realize indirectly by means of a high-speed infrared cam- 
era. Variations in the temperature over time are subsequently 
obtained as the mean value over the pixel intensity of every 
infrared snapshot. We start unraveling the dynamics of the ob- 
served laser droplet generation via the embedding theorem, 23 
which enables the reconstruction of the phase space from a 
single observed variable, thereby laying foundations for fur- 
ther analysis. We use the mutual information 24 and the false 
nearest neighbor method 25 ' 26 to obtain optimal embedding pa- 
rameters for the phase space reconstruction. Subsequently, 
we apply a determinism 18 and a stationarity 21 test to verify 
that the observed behavior is indeed a consequence of deter- 
ministic dynamics and that all the parameters were held con- 
stant during data acquisition. After establishing that the stud- 
ied temperature recording originates from a deterministic and 




FIG. 1: Schematic presentation of the experimental setup. The tem- 
perature is measured indirectly by means of a high-speed infrared 
(IR) camera (see main text for details). 



stationary laser droplet generation, we calculate the spectra 
of Lyapunov exponents, 27 ' 28 whereby a positive largest Lya- 
punov exponent— and a negative divergence both point to- 
wards the fact that the observed dynamics is deterministic ally 
chaotic with an attractor as solution in the phase space. We 
also outline potential industrial applications of the process and 
give pointers towards its appropriate mathematical modeling. 
It is notable that the number of droplet-based technologies has 
increased substantially in recent years. 30 From this point of 
view laser droplet generation has an industrial potential es- 
pecially in joiningj^L^ 2 . where accurate mathematical models 
would be paramount for further application developments. 

The paper is structured as follows. Section II is devoted 
to the accurate description of the experimental setup and the 
acquisition of the time series. Section III features results of 
nonlinear time series analysis, while in the last Section we 
summarize the paper and outline potential implications of our 
findings. 

II. EXPERIMENTAL SETUP 

The laser droplet generation phenomenologically consists 
of two phases. In the first phase a laser pulse is used as a 
source of energy to melt the end of a vertically placed metal 
wire. From the molten end a pendant droplet is formed due 
to the action of surface tension and gravity force. Because the 
surface tension drags the pendant droplet up the wire the lat- 
ter has to be properly fed. The second phase encompasses the 
detachment of the pendant droplet from the tip of the wire. To 
achieve this the surface tension force needs to be overcome, 
which in our case happens as a result of droplet mass growth. 
This scenario can be referred to as spontaneous dripping. No- 
tably, a metal droplet can be used in different manufactur- 
ing applications. The most promising one is droplet joining, 
where a molten droplet is placed onto the joining spot^I^ 
The heat content of a droplet is sufficient to produce a high- 
temperature weld, whereas the volume of the droplet can be 
used to fill gaps or bridge dimensional tolerances. Other po- 
tential applications include the generation of 3D structures 
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FIG. 2: A sequence of snapshots from the infrared camera showing 
growth of the pendant droplet and its detachment. The time increases 
from the upper left to the lower right panel. We also provide a sup- 
plementary video file from the infrared camera, showing the process 
of laser droplet generation in continuous time (see snapmovie.avi). 



by means of a selective deposition of droplets into layers as 
well as micro casting. The most common and important pro- 
cess underlying these technologies is laser droplet generation. 
However, for an effective optimization and control of the pro- 
cess it is essential to know its dynamics. We aim to determine 
this from experimental data. 

In order to study the dynamics of laser droplet generation 
we have developed an experimental system that is schemati- 
cally depicted in Fig. Q] The main parts of the experimental 
setup are the Nd:YAG pulse laser, the opto-mechanical ele- 
ments, the wire feeder, and the infrared camera. The Nd:YAG 
laser is used for generating laser pulses with a wavelength of 
1 .06/im. The maximal laser pulse power is 8kW and the pulse 
duration is between 0.3ms and 20ms. The maximal pulse rep- 
etition rate is 300Hz with an average power of 0.25kW. To 
assure uniform heating of the wire and process symmetry, the 
laser light is divided into three equal laser beams. By means 
of the opto-mechanical elements the beams are distributed 
equiangular along the wire circumference and perpendicularly 
focused onto the wire's surface. The wire is vertically fed via 
a controlled wire feeder. Maximal acceleration and velocity of 
the wire are 20m/s 2 and 0.3m/s, respectively. The wire con- 
troller is also applied to synchronize the triggering of the laser 
pulses with the stepwise wire feeding. Since the temperature 
is the most important variable of the process it was indirectly 
measured by means of a high-speed infrared camera. Given 
the properties of the light emitted at the wire end and the pen- 
dant droplets, the snapshots were acquired at wavelengths be- 
tween 3.5/im and 5/im. 

Based on the description of the experimental setup it is ob- 
vious that there are a number of parameters that can influ- 
ence the process of laser droplet generation. For a selected 
wire material, however, the most important ones are the laser 
pulse and wire feeding parameters. Proper estimates for the 
laser pulse parameters can be obtained based on the analysis 
of the heat balance of a molten pendant droplet. 32 Parameters 
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FIG. 3: The time series, capturing the process of laser droplet gen- 
eration via the pixel intensity of high-speed infrared snapshots (see 
Fig. [2}- The main panel shows the time series Xi rescaled to the unit 
interval with an integer time scale. The inset shows the series before 
(red line) and after (black line) Wiener filtering, removing the high- 
frequency noisy component that is due to the infrared imaging. No 
additional noise filtering has been made prior to further analysis. 

of the wire feeding generally depend on the dynamics of the 
droplet. However, an average wire feed velocity can be esti- 
mated based on the desired droplet volume. 32 Here we have 
used a nickel wire having diameter 0.6mm. Other parameters 
and setup details are as follows: A rectangular laser pulse with 
power 1440W, duration 12ms and frequency 3Hz was used. 
During the laser pulse, the wire was fed by a triangular veloc- 
ity profile with a maximal velocity of 0.3m/s. The motivation 
behind the usage of a triangular profile of the feeding velocity 
is two-fold. First, it is simple and transparent enough to be 
easily implement experimentally, and second, due to its sim- 
plicity, it is as non-invasive on the inherent dynamics of the 
droplet generation as possible, in particular, allowing swift 
adjustments in accordance with the droplet growth and sub- 
sequent detachment. The sampling frequency of the infrared 
camera was 1428Hz at snapshot size of U x V — 32 x 64 pix- 
els. A short sequence of snapshots is shown in Fig. |2] where 
an example of droplet growth and subsequent detachment is 
depicted (see also the supplementary video snapmovie.avi). 
Finally, the spatiotemporal temperature field was converted 
into a single scalar time series by calculating the mean value 
of the pixel intensity of every snapshot according to: 

1 U V 

U— 1 V— 1 

In Eq. (fl3 x\ ,v is the z-th snapshot value of the (u,v) pixel, 
and thus is as a proxy of the local temperature T(u, v, t) = 
x"' v , where t = idt [see also Eq. (0]. The resulting time 
series is, rescaled to the unit interval and de-noised by means 
of a Wiener filter, shown in Fig. [3] Upon visual inspection of 
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FIG. 4: Determination of the proper embedding delay via the mutual 
information method. The first minimum occurs at r ~ 250 (blue 
dashed line), which we will use in all subsequent calculations. 



the time series, lower and higher frequency oscillations can be 
inferred, which can be linked nicely with the two-phase pro- 
cess of laser droplet generation. Namely, the lower frequency 
oscillations correspond to droplet volume (mass) and temper- 
ature growth, which is followed by a sudden drop of the signal 
amplitude due to the droplet detachment. 



III. TIME SERIES ANALYSIS 

We start the time series analysis by applying the embedding 
theorem, 23 which states that for a large enough embedding 
dimension m the delay vectors 

z(l) = [Xi, X{+ T , X{+2t, ■ ■ ■ i x i+(m-l)T] (2) 

yield a phase space that has exactly the same properties as 
the one formed by the original variables of the system. In 
Eq. (O variables Xi, Xi +T , xi + 2 T ,- ■ ■ , ^i+( m -i)T denote val- 
ues (rescaled to the unit interval for simplicity) of the indi- 
rectly measured temperature at times t — idt, t = (i + r)dt, 
t = (i+2r)d<,. . . ,t = [i+(m— l)r]di, respectively, whereby 
r is the embedding delay and dt is the sampling time of data 
points equaling 7- 10~ 4 s. Altogether, the examined time series 
consists of i = 1,2,..., 66732 data points. 

While the implementation of Eq. (O is straightforward, we 
first have to determine proper values for the embedding pa- 
rameters m and r. For this purpose, the mutual information 24 
and the false nearest neighbor method 25 can be used, respec- 
tively. Since the mutual information between Xi and Xi+ T 
quantifies the amount of information we have about the state 
Xi+r presuming we know Xi,^- Fraser and Swinney — pro- 
posed to use the first minimum of the mutual information 
as the optimal embedding delay. Results presented in Fig. [4] 
show that the mutual information I(t) has the first minimum 
at t pa 250. The false nearest neighbor method, on the other 
hand, relies on the assumption that the phase space of a de- 



FIG. 5: Determination of the minimally required embedding di- 
mension. The fraction of false nearest neighbors (fnn) drops close 
(< 0.01) to zero (blue dashed line) at m = 5, which we will use in 
all subsequent calculations. 



terministic system folds and unfolds smoothly with no sud- 
den irregularities appearing in its structure. By exploiting this 
assumption one comes to the conclusion that points that are 
close in the reconstructed embedding space have to stay suf- 
ficiently close also during forward iteration. If a phase space 
point has a close neighbor that does not fulfil this criterion 
it is marked as having a false nearest neighbor. As soon as 
to is chosen sufficiently large, the projection effects due to 
a mapping of the time series onto a space with too few de- 
grees of freedom should disappear, and with them the fraction 
of points that have a false nearest neighbor (fnn) should con- 
verge to zero. 25 Note that the method implicitly assumes that 
a deterministic time series is given as input. This, however, 
cannot be taken for granted, and indeed a simple extension of 
the originally proposed false nearest neighbor method 2 ^ can 
be used also as a determinism test. Here we employ the clas- 
sical algorithm proposed by Kennel et al— and use the deter- 
minism test due to Kaplan and Glass 18 . Results of the false 
nearest method are presented in Fig. [5] showing that fnn — > 
at m = 5. We will thus use r = 250 and m = 5 as input for 
Eq. d2J m what follows. 

Having all the parameters at hand for reconstructing the 
phase space from the observed variable (see left panel of 
Fig. [6j, we can proceed by employing the determinism test 
proposed by Kaplan and Glass. 18 The test is simple but effec- 
tive, measuring average directional vectors in a coarse-grained 
embedding space. The idea is that neighboring trajectories in 
a small portion of the embedding space should all point in the 
same direction, thus assuring uniqueness of solutions in the 
phase space, which is the hallmark of determinism. To per- 
form the test, the embedding space has to be coarse grained 
into equally sized boxes. The average directional vector per- 
taining to a particular box is then obtained as follows. Each 
pass p of the trajectory through the fc-th box generates a unit 
vector e p , whose direction is determined by the phase space 
point where the trajectory first enters the box and the phase 
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FIG. 6: Determinism test. Left panel features the reconstructed phase 
space using r = 250 and m — 5, while the right panel shows the 
pertaining approximated directional vector field. Determinism factor 
of the phase space according to Kaplan and Glass 18 is k = 0.9. 
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FIG. 7: Stationarity test. The whole time series was partitioned 
into 26 non-overlapping segments each occupying 2500 data points. 
The color map displays average cross-prediction errors S g h in depen- 
dence on different segment combinations. 



space point where the trajectory leaves the box. The average 
directional vector Vfe through the k-th box is then 

n 

V k = n- 1 J2 e P ( 3 ) 

P =i 

where n is the number of all passes through the fc-th box. 
Completing this task for all occupied boxes gives us a di- 
rectional approximation for the vector field. If the time se- 
ries originates from a deterministic system, and the coarse 
grained partitioning is fine enough, the obtained directional 
vector field should consist solely of vectors that have unit 
length. Hence, if the system is deterministic, the average 
length of all the directional vectors k will be close to one. 
The determinism factor pertaining to the five-dimensional em- 
bedding space presented in Fig. [6] that was coarse grained 
into a 12 x 12 x ... x 12 grid is n = 0.9, which con- 
firms the deterministic nature of the studied time series. The 
two-dimensional projection of the directional vector field is 
shown in the right panel of Fig. [6] It is also informing to 
generate surrogates 34 from the studied series, for example by 
employing the iterative procedure proposed by Schreiber and 
Schmitz, 35 albeit with the cautionary note that surrogates en- 
able only the rejection (or acceptance) of a given null hy- 
pothesis. As such, they do not allow more far-reaching con- 
clusions on the role of stochasticity in the examined series, 
and thus cannot be used as a substitute for the preceding de- 
terminism test. In our case the determinism factor drops to 
0.48 < k < 0.61 (based on 20 generated surrogates), which 
thus rejects the null hypothesis that the laser droplet genera- 
tion is a stationary Gaussian linear process that has been dis- 
torted by a monotonic, instantaneous, time-independent non- 
linear function. 

It remains of interest to verify if the laser droplet gener- 
ation is a stationary process. To this purpose, we apply the 



stationarity test proposed by Schreiber.- 1 In general, station- 
arity violations manifest so that various non-overlapping seg- 
ments of the time series have different dynamical properties. 
Since linear statistics, such as the mean or standard data de- 
viation, usually do not posses enough discrimination power 
when analyzing irregular signals, nonlinear statistics have to 
be applied. One of the most effective has proven to be the 
cross-prediction error statistic. The idea is to split the time 
series into several short non-overlapping segments, and then 
use a particular data segment to make predictions in another 
data segment. By calculating the average prediction error 
S g h when considering points in segment g to make predic- 
tions in segment h, we obtain a very sensitive statistics capa- 
ble of detecting minute changes in dynamics, and thus a very 
powerful probe for stationarity. 21 Results presented in Fig. [7] 
were obtained by dividing the whole time series into 26 non- 
overlapping segments of 2500 points, thus yielding 26 2 com- 
binations to evaluate Sgh. Since the cross-prediction errors 
are uniformly spread across the whole g — h plane, i.e., none 
of the segments is an exceptionally bad (or good) source of 
data to make predictions in the other segments, we can re- 
fute non-stationarity in the proposed laser droplet generation 
experiment. An interesting feature of the stationarity test pre- 
sented in Fig.|7]is also the emergence of white diagonals, ap- 
pearing parallel to the main white diagonal. While the latter 
is expected because there g — h (the segment used for mak- 
ing predictions is also the one we test them against), the other 
white diagonals appear due to the pseudo-periodicity, which 
can be inferred from the outlay of the time series (see Fig. [3}. 
In particular, although the growth and subsequent detachment 
of droplets is obviously not a strictly periodic process, the sim- 
ilarity of the individual periods of droplet growth, as well as 
the similar lengths of time spans between consecutive droplet 
detachments, nevertheless lead to the recurrent emergence of 
white diagonals also beyond g = h. This, in turn, can also be 
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FIG. 8: Spectra of Lyapunov exponents determined using radial basis 
functions for the approximation of the flow. From top to bottom the 
lines depict the convergence of the largest (Ai) to the smallest (As; 
most negative) Lyapunov exponent as a function of the discrete time 
i. The lowest line (pink dashed) shows the sum of all five exponents, 
i.e., the divergence A = 2j=i m ^i- ^ nnear fit towards the end 
of the curves gives Ai = (3.2 ± 0.1)s _1 , A 2 = (0.0 ± 0.1)s" 1 and 
A = —(145 ± 3)s _1 . Note that the vertical axis has a break. 

seen as evidence supporting the determinism in the examined 
time series. 

Finally, we calculate the spectra of Lyapunov exponents \j 
where j — 1,2, ... ,m, knowing with reasonable certainty 
that the obtained results are due to deterministic nonlinear dy- 
namics rather than noise or varying systems parameters during 
data acquisition. We employ radial basis functions for the ap- 
proximation of the flow in the phase space. Using the phase 
space reconstruction parameters obtained above, M = 10 
nearest neighbors of each z(i) to make the fit, and the stiff- 
ness parameter r = 7, 27 the exponents change their sign upon 
time reversal of the flow and converge robustly as the number 
of iterations increases. Figure [8] features the individual con- 
vergence curves, from which we obtain Ai = (3.2 ± 0.1)s _1 , 
A2 = (0.0±0.1)s _1 , and the divergence as the sum over all Xj 
equal to A = —(145 ± 3)s~ 1 . From the positive largest Lya- 
punov exponent, the vanishing second Lyapunov exponent, 
and the negative divergence, we can conclude that the dynam- 
ics of laser droplet generation is deterministically chaotic, and 



that there exists a stable attractor in the phase space to which 
any given cloud of initial condition converges in time. 



IV. SUMMARY 

We have proposed an experimental setup with the aim of 
determining the dynamics of laser droplet generation. Using a 
high-speed infrared camera, we have indirectly measured the 
spatiotemporal profile of temperature around the molten end 
of the wire and the pending droplets. Subsequently, the time 
series was obtained as the mean value over the pixel intensity 
of every infrared snapshot, and analyzed systematically with 
methods of nonlinear time series analysis. After reconstruct- 
ing the phase space from the observed variable, we have veri- 
fied that the later has properties that are typical for determin- 
istic and dynamically stationary systems. We have shown that 
the minimally required embedding dimension is five, which 
altogether suggests that it would be justified to mathematically 
model the process of laser droplet generation with no more 
than five first-order ordinary differential equations. Also, we 
have determined the whole spectra of Lyapunov exponents by 
approximating the flow in the phase space with radial basis 
functions. Our calculations revealed that the largest Lyapunov 
exponent is positive, the second is zero, while the divergence 
is negative, thus obtaining strong indicators that the observed 
dynamics is deterministically chaotic with an attractor as so- 
lution in the phase space. Thus, although the laser droplet 
generation is governed by additional physical phenomena, in- 
cluding light-metal interaction, heating and phase transitions, 
the dynamics of the process is similar to the one observed in 
traditional dripping faucet experiments. In addition, the pre- 
sented results indicate that nonlinearity is an innate ingredient 
of laser droplet generation, which should be taken into ac- 
count in future modeling and controlling attempts. We hope 
that the study will be of value when striving towards a deeper 
understanding of the examined process and its integration into 
outlined industrial applications. 
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